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Abstract 

We examine the stability of wireless networks whose users are distributed over a 
compact space. Users arrive at spatially uniform locations with intensity A and each 
user has a random number of packets to transmit with mean 0. In each time slot, an 
admissible subset of users is selected uniformly at random to transmit one packet. A 
subset of users is called admissible when their simultaneous activity obeys the prevailing 
interference constraints. We consider a wide class of interference constraints, including 
the SINR model and the protocol model. Denote by \i the maximum number of users 
in an admissible subset for the model under consideration. We will show that the 
necessary condition A/3 < (i is also sufficient for random admissible-set scheduling to 
achieve stability. Thus random admissible-set scheduling achieves stability, if feasible to 
do so at all, for a broad class of interference scenarios. The proof relies on a description 
of the system as a measure- valued process and the identification of a Lyapunov function. 



1 Introduction 

The present paper examines the stability of a broad class of wireless networks whose users 
arrive to a compact space H. Time is slotted, and users arrive at H according to some spatial 
stochastic process with mean A per time slot. Users independently take their locations in H 
at random according to the uniform distribution. Each user has a random number of packets 
to transmit, generally distributed with mean (3, and can transmit at most one packet per 
time slot. 

In each time slot, we select a set of users for transmission from all admissible sets uniformly 
at random. A subset of users is called admissible when their simultaneous activity obeys 
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the prevailing interference constraints. In practice, the relevant interference constraints 
depend on various system-specific properties, such as the propagation environment and the 
operation of the physical and medium access layers of the network. In the present paper we 
therefore adopt generic feasibility criteria, which in particular cover both the SINR (Signal- 
to- Interference-and-Noise Ratio) model and the protocol model as two canonical models for 
interference. 

Let \x be the maximum number of users in an admissible set. It is clear that A/3 < /i is 
a necessary condition for stability: The mean number of packets that arrive per time slot 
should be no larger than the maximum number of packets that can be transmitted simul- 
taneously. The main result we show is that this necessary condition is also nearly sufficient 
for stability. Specifically, the Markov chain describing the evolution of the system is then 
positive Harris recurrent, implying that the network will be empty infinitely often. Thus 
random admissible-set scheduling is a highly robust strategy in that it achieves stability, if 
feasible to do so at all, for a wide range of interference scenarios. 

For wireless networks we are primarily interested in interference constraints that become 
looser when users are further apart and spaces such as the one- and two-dimensional torus 
and the ordinary sphere. Our results, however, hold in more generality, which is why we 
present our results in terms of a particle system on a compact space H, where particles (or 
packets) arrive to the system according to a spatial stochastic process with rate A and in 
batches with mean /3. The arrival location is independent of other users and uniform on H. 
As described above, we investigate stability in the context of a model that combines a 
scheduling discipline operating under interference constraints and a continuous spatial set- 
ting. While these two elements have each been considered in isolation before, the present 
paper is, to the best of our knowledge, the first to capture both features in conjunction. In- 
deed, stability of wireless networks has been widely studied in the literature, see for instance 
Bonald and Feuillet [2], Bordenave et al. |3j and Wu et al. |12j . These papers restrict the 
attention though to discrete topologies and interference constraints such that the system 
can be represented as a conflict graph. Our model does not allow such a representation due 
to the continuum of locations, and hence these results are not directly applicable to our 
problem. Stability of queueing networks in continuous space is investigated in Altman and 
Levy pQ, Robert |T0] and Leskela and Unger [7J. These papers prove stability of networks in 
which only one user is allowed to transmit at a time. In contrast, the present paper focuses 
on the more complex situation of simultaneous transmissions as governed by a scheduling 
discipline. 

While a discrete network structure is a reasonable assumption in case of a relatively small 
number of long-lived sources, it is less suitable in case of a relatively large number of 
short-lived flows. The latter scenario is increasingly relevant as emerging wireless networks 
support traffic generated by massive numbers of nodes which each individually may only 
engage in sporadic transmission activity. The continuous spatial setting also provides useful 
insights in the scaling behavior of discrete topologies as the number of nodes grows large. 
From a methodological perspective, the continuous spatial setting involves major additional 
challenges compared to a discrete network structure. Since users reside in a continuum of 
locations, the evolution of the system cannot be represented in terms of a Markov chain with 
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a finite state space, and we therefore introduce a measure-valued process as a description 
of the system. In order to prove stability, we identify a Lyapunov function which has a 
negative drift for all but a 'small' set of states, so that the Markov chain is positive Harris 
recurrent. 

The remainder of the paper is organized as follows. In Section [2] we show how the evolution 
of the system may be described in terms of a Markov chain with a measure-valued state 
space. The main stability result is presented in Section [3], along with an interpretation and 
discussion of its ramifications. In Section [4] we provide the proof for our main result, and 
in particular identify a Lyapunov function which has negative drift for all but a 'small' set 
of states, and plays a critical role in the proof. In Appendix [A] we recall various useful 
definitions and collect some preliminaries that are needed in order to apply the Foster- 
Lyapunov approach for our specific Markov chain. 

2 Model description 

Consider a compact space H C M n for some n < oo. Denote by A(t, B) the number of 
particles arriving during the i-th time slot in B C H. Particles arrive in batches, and the 
batch size has a general non-negative discrete distribution with mean (5, independent of the 
sizes and locations of other batches. We assume batches have a size of at most one with 
positive probability. Batches arrive at locations uniformly distributed on H, independent 
of the locations of other batches. The number of batches that arrive during a time slot has 
a general non-negative discrete distribution with mean A and is independent of the number 
of batches that arrive in other time slots. That is, the numbers A(t,H), t = 1,2, ... , are 
i.i.d. copies of a non-negative random variable A with E{yl} = A/3. Further we assume 
that E{Alog(A)|A > 0} < oo. Note that, because batches arrive uniformly, the expected 
number of particles to arrive to a subspace B C H in one time slot is given by E {^4} v{B), 
where v is some measure on H such that, without loss of generality, v{H) = 1, representing 
for example the volume or surface area. 

We denote the number of particles in the space H at the start of the i-th time slot by Y(t), 
with Y(t) = (Y(t, B), B C H) and Y(t, B) denoting the number of particles residing in the 
subspace B at the start of the t-th time slot. The state space of this process is denoted 
by ^ and consists of all finite counting measures on H. So, when y G y(B) denotes the 
number of particles residing in B. In particular, y({x}) denotes the number of particles at 
x G H. 

At the start of every time slot an admissible subset of particles will be removed. Here, 
z G ^ is called a subset of y G ^ if z({x}) < y({x}),Vx G H. To decide whether a set is 
admissible we define a function F : ^ — > {0, 1} having the following properties 

1. F(yo) = 1, where yo denotes the empty configuration, yo(H) = 0. 

2. F(y) = 1 only if y({x}) G {0, 1}, Vx G H. 

3. There exists a partition of H, i.e. disjoint subspaces P{, i = 1, . . . , K, for which v{Pi) = 
i,Vi G {1,...,K} and P% = H, such that F(y) = 1 with y({x}) = 1 for 
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some x G Pj implies that y(Pi \ {x}) = 0. Furthermore, taking y G ^ such that 
y(P m+ j) = l,Vj G and y(.ff \ Uj=i ^W+j) = 0, with m < K such that 

m mod ^ = 0, implies that F{y) = 1. Here [i = max{y(H) : F(y) = 1}, i.e. the 
maximum number of particles in an admissible set, and K is a multiple of \x. 

We call !/ 6 $ admissible if and only if F(y) = 1. 

Property [T] ensures that the empty set is admissible, such that there always is at least one 
set that is admissible. Property [2] states that at most one particle of the particles located 
at some location x G H can be removed in a time slot. Property [3] states that the space H 
can be partitioned using a finite number of equally-sized disjoint sets such that at most 
one particle in each set can be removed in a time slot. This property further states that 
this partitioning is such that taking any combination of elements of maximal size (^) from 
certain sets will always give an admissible set, this is to avoid having almost surely no 
admissible subset of particles of maximum size at any point in time. Finally, Property [3] 
states that the partitioning is such that there are K/ fi disjoint sets of sets with this property, 
i.e. the property makes sure that taking any combination of elements of maximal size from 
these sets will always give an admissible set. Note that fi is defined implicitly. 
We will now give an example of a function which satisfies the above properties. The model 
is this example is the so-called protocol model. 

Example 2.1. Consider a unit circle, use the interval [0, 1) to denote points on the circle 
and let v be the Lebesgue measure. Particles can be removed simultaneously whenever the 
distance between these particles is at least r, i.e. F(y) = 1 if y({x}) G {0, 1}, Vx G [0, 1) and 
y({x})y({w}) > only if D(x, w) > r Vx ^ w G [0, 1), with D(x, y) = min(|x— y\, 1 — \x— y\). 
Properties [T] and [2] immediately follow from this description. To verify property [3] take 
K > 2fi/ (1 — /ur) and K a multiple of fi = |l/ r J • Now, for i = 1, . . . , K, take 



We then see that with these sets property [3] is satisfied whenever 1/r is non-integer. For 
property [3] to hold if 1/r is integer- valued we do not allow sets of size 1/r to be removed. 
Note that the probability that a set of size 1 jr is removed is almost surely zero in any time 
slot if 1/r is integer- valued. We then see that property [3] is satisfied with the above sets Pi 
and fi = [1/r] — 1. 

Let x{y) be the set of all subsets of y, i.e. 



and let R(t, Y(t), B) be the number of particles removed from B C H in the t-th time slot, 
given the configuration, Y(t). An admissible subset of particles is selected uniformly at 
random. Hence, given Y(t) = y, R(t,y) = z with probability 




X (y) = ^f: z{{x}) < y({x}),Vx G H} 



f(z) n 



x<=H,z({x})>0 
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where R(t,y) = (R(t,y, B), B C H). Note that for z x(u) we have F(z) = or 3x € H 
such that z({x}) = 1 and y({x}) = 0, so this probability is always zero in this case. Further 
note that, here and in the remainder of this paper, the value of the empty product is defined 
as 1, as is usually done in the literature. 
The evolution of Y(t, B) is then described by the recursion 

Y{t, B) = Y(t — 1,B) + A(t -1,B)- R(t, B). 

Further, R(t, B) depends on the number of particles just before the start of the i-th time 
slot, so 

Y{t~, B) = Y(t -l,B) + A(t -l,B). 

From this description it follows that (Y(t))t£N is a Markov chain. We will equip the state 
space of the Markov chain, vf, with the smallest cr-field with respect to which the 

map y — > y{B) is measurable for any Borel set B C H. That is, we equip ^ with the Borel 



cr-field as we will prove in Lemma A. 2 



3 Main result 

The following theorem states the main result of this paper and is a shorter version of 
Theorem 4.5 which is proven at the end of Section [4} 



Theorem 3.1. Assume A/3 < fi. Then, the Markov chain (Y(t))t(z?>} is positive Harris 
recurrent. 

This theorem states that, starting from the empty configuration, Y(t, H) = for infinitely 
many values of t, i.e. the Markov chain will be in the empty configuration infinitely often. 
We thus see that random admissible-set scheduling achieves maximal stability, independent 
of the specific function F. 



The result of Theorem 3.1 may be interpreted as follows. Suppose that the total number 
of particles in the system is large. Then there will be a large number of admissible sets 
of size fi, assuming that the particles are sufficiently dispersed across the network and not 
concentrated in a few dense areas. In fact, the number of admissible sets of size /i will 
be overwhelmingly large compared to the number of admissible sets of smaller size. By 
virtue of random-admissible set scheduling, one of the admissible sets of size fi will then be 
selected with high probability. Thus, the expected number of removed particles will exceed 
the expected number of arriving particles, provided A/3 < /j,, implying a reduction in the 
expected number of particles in the system, and preventing the number of particles from 
growing without bound. 

As the above heuristic explanation indicates, it is crucial for the particles to be sufficiently 
spread out and not be clustered in a few hot spots. In order to obtain a rigorous proof, it 
will hence not suffice to just consider the total number of particles, but in fact be necessary 
to keep track of their individual locations. 
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Figure 1: Random- admissible set scheduling and maximal scheduling with priorities for 
( = 0.5, r = 0.49 and A = 1.95. 



It is worth emphasizing that we consider a particle-based version of random- admissible 
set scheduling rather than a node-based incarnation, in the sense that the strategy selects 
among sets of particles rather than sets of nodes. While this distinction is immaterial when 
the measure u(-) is an absolutely continuous density, the issue does become relevant when 
the measure z/(-) has mass in discrete points. In the latter case, it may readily be concluded 
that a node-based version of random admissible-set scheduling may fail to achieve maximum 
stability. Interestingly, this observation contrasts with the fact that the node-based version 
of the celebrated MaxWeight scheduling strategy guarantees maximum stability when the 
measure u(-) is a purely discrete distribution, whereas a particle-based version may fail to 
do so p]. 

Also, while the spatial dispersion of particles under random admissible-set scheduling is 
intuitively plausible, it is certainly not obvious. This is perhaps best illustrated by the fact 
that the result of Theorem 3.1 may not necessarily hold for seemingly similar but subtly 
different scheduling disciplines. 

As an example consider the situation of Example 2.1 Instead of selecting an admissible set 
at random, now consider the scheduling discipline that gives priority to particles that are 
closest, in anticlockwise direction, to a certain point £ on the circle. Obeying this priority 
rule, we select as many particles as possible to remove in a certain time slot. That is, the 
particle closest to the given point gets removed, the particle closest to the point and at least 
a distance r away from the first particle gets removed, and so on until no particle can be 
selected anymore. We call this service discipline maximal scheduling with priorities. 
Figure [T] shows a simulation result for both scheduling disciplines with Q = 0.5, r = 0.49 and 
A = 1.95 starting from an empty configuration and running for 10 6 time slots. That is, the 
figure gives a realization of Y(t, H) given that 1^(0, H) = for t = 1, . . . , 10 6 . We see that at 
the start of the simulation the number of particles in the system with random-admissible set 
scheduling grows faster than the number of particles in the system with maximal scheduling 
with priorities. This is because maximal scheduling with priorities always selects a subset of 
maximum size obeying the priority rules, whereas random- admissible set scheduling always 
selects admissible sets of a small size with a certain probability, which gets lower as the 
number of particles in the network grows. More importantly, after some time the number of 
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Figure 2: Terminal configuration of the simulation of Figure [TJ 

particles in the system with random- admissible set scheduling settles around an equilibrium 
value whereas the number of particles in the system with maximal scheduling with priorities 
keeps on growing linearly. This suggests that maximal scheduling with priorities is not stable 
while random-admissible set scheduling is stable for the chosen parameters. The latter will 
be proven in the next section. 

Figure [2] shows the terminal configuration of the simulation of Figure [TJ i.e. it gives a 
realization of Y(W 6 ,B) given that Y(0, 1) = for B = [0, s), < s < 1, for both scheduling 
disciplines. For random admissible-set scheduling we see that the number of particles in the 
interval [0, s) is roughly linear in s, indicating that the particles are evenly spread out over 
the circle. For maximal scheduling with priorities we observe that the number of particles 
in [0, s) slowly increases with s up to approximately s = 0.48, after which the number of 
particles in the system steeply rises up to s = 0.5. For s > 0.5 = C, the number of particles in 
the interval [0, s) is (almost) constant, implying that virtually no particles are located in the 
interval [0.5,1). Note that particles in the interval [0.48,0.5) have the lowest priority and 
hence are, whenever they are allowed to, almost always removed simultaneously with other 
particles, as there are quite some particles in the system and outside this interval. However, 
the particles that are allowed to be removed simultaneously with particles in [0.48,0.5) are 
also allowed to be removed simultaneously with some particles in [0.5,0.52), who have the 
highest priority. So we infer that the particles in this system are clustered in [0.48,0.5), 
and that too high a fraction of the time (larger then 0.02A) no particle in this interval is 
removed, making the system unstable. 



4 Proof 

In this section we provide the proof of Theorem 3.1 As mentioned earlier, the proof relies on 



the Foster-Lyapunov criteria and involves the identification of a function which has negative 
drift for all but a small set of states. Appendix [A] contains several useful definitions and 
preliminaries that are needed to apply the Foster-Lyapunov approach for our specific Markov 
chain. 

In order to define the Lyapunov function, we use a partitioning Pi, ... , Pk of the space H 
for which Property [3] of the function F holds. Let Xk(y) be the number of particles residing 
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in the k-th region, P^, given configuration y G X I / , for k G fC, where /C = {1, . . . , K}. 
Let f2 = V(IC) be the collection of all subsets of {1,. . . , K} and let Q(y) be the subsets 
containing particles, i.e. Q(y) = {S G ft : Xk(y) > l,Vfc G S}. A subset S G is called 
'guaranteed' if any subset of particles, with exactly one residing in each of the regions 
contained in S, is admissible, regardless of the exact locations within each of the regions. 
Let C S] be the collection of all 'guaranteed' subsets and, again, let 6(y) be the subsets 
containing particles, 0(y) = {S G G : Xfc(y) > 1, Vfc G S}. Further, denote by qs(y) the 
probability that the particles that get removed belong to the subset of regions S G ft, given 
that the system is in state y, with Pk{y) = ^SGtt-SBkQsiv) the probability that a particle 
in the k-th region gets removed given configuration y G ^f. 
For S £ Q, denote 

w s(y) = J| Xk{y). 

keS:x k (y)>l 

Further define 

B(e) = jy G ^ : w{y) > 

where e > and 

w{y) = max ws(y)- 
see( y ) 

In the next lemma we will show that the value of Ylkes ^°s( x k), for the set of particles S 
selected with random admissible-set scheduling, is close to the maximum possible value over 
all admissible sets with high probability for all states y G B(e). 

Lemma 4.1. For all states y G B{e) we have, 
qs ( y "> lo &( w s(y)) > (! - e ) log(w(y)). 

Proof. The proof proceeds along similar lines as in [2J, [3]. 
Define 

T(y) = {S G Q(y) : log(^(y)) > (1 - |) log(«;(y))}. 

Then, 

<ls(y)log(ws(y)) > (1 - |)log(iu(y)) ^ gs(y). (1) 
sen(y) ~ SeT(y) 

For any 5 € O, let i>s(y) be the number of admissible subsets of particles of size |5| with 
exactly one residing in each of the regions contained in S, with the convention that v$(y) = 1 
for all y G Because every admissible subset can have at most one particle residing in 
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each region, there is exactly one S G f2 for which this subset is counted in vs(y)- Thus the 
total number of admissible subsets of particles is given by X^Tefi v t(u) an d, as an admissible 
subset of particles is selected uniformly at random, 

, \ v s (y) v s (y) 

qs(y) - 



Further observe that vs(y) < ws{y) for all S G Q, with equality for all S G 6(y). 
Thus 

s ^ {y) Esen(y) v s(y) Esee(y) w s(y) w{y) 



The latter quantity is less than | for all states y G B(e), and thus 

e 
2 



2 

E «s(y)>i-|- ( 2 ) 

Combining the lower bounds ([!]) and ([2]), we obtain 

E qs ^ l °g( w s(y)) > (1 - |) log(w(i/))(l - |) > (1 - e) logMy)) 

for all states y G B(e). □ 

For i = 1, . . . , K define the set Si G Q by Si = { |~^] [i — j, Vj G {0, . . . , fj, — 1}. Note that, 
by definition of F, Si G 0. Further define Si(y) = {j G Si(y) : Xj(y) > 1} and note that 
Si(y) G G(y). 

Also define the nonnegative function V : * — > M. by 

= E x k(y) l og(x k (y)). 

and the function G : f 1 by 

G(y)= E log(^(y))[E{^ fc }- Pfc (y) 
fe:a!fc(i/)>l 

where A k denotes the number of arrivals in the k-th region, so E{j4/%} = 
Note that 

K 

v (y) = ^2x k (y)log(max(x k (y),l)) 



k=l 



and 



V(y)= E x k(y)^og(x k (y)). 

k-x k {y)>2 
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Further observe that the function V(-) only depends on y through the values of x k (y). 
However, the Xfc(y)'s do not constitute a Markov chain, and hence we need to treat V(-) 
as a function of the full state description y in order for the Foster-Lyapunov approach to 
apply. 

We now first find the relation between the drift of V(y) and G(y), where the drift of V(y) 
is defined by 

AV(y) =E{(Y(t + l))\Y(t) = y} - V(y). 

After that we will find an upper bound for G(y). 

Lemma 4.2. AV(y) = G(y) + G^tj/), with G^y) a bounded function. 

Proof. Remember that at most one particle can be removed from each region in every time 
slot. We thus get 

E{V(Y(t + l))\Y(t)=y} 

K 

= Elj2 x k(Y(t + 1)) log(max(x fe (Y(i + 1)), 1))| Y(t) = y) 
k=l 

K 

= ^E{x fe (y(t + 1)) log(max(x fc (Y(t + 1)), l))\r(t) = y) 
k=i 

Y Pk{y)^{{x k {y) + A k (t)-l)\og(x k (y) + A k (t)-l)} 

k:x k (y)>2 

+ Y (l-Pk(yM{(x k (y)+A k {t))log(x k (y)+A k {t))} 

k:x k {y)>2 

+ Y Pk(y)W{A k (t)>l}E{(l + A k (t)-l)log(l + A k (t)-l)\A k (t)>l} 

k:x k (y)=l 

+ Y (l-Pfc(j/))E{(l + A fc (t))log(l + A fc (*))} 

k:x k (y)=l 

+ Y ^{A k (t)>l}E{A k (t)log(A k (t))\A k (t)>l}, 

k:x k (y)=0 
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We further have 

Pk(y)^{(x k (y) + A fc (t) - 1) log(x k (y) + A k (t) - 1)} 

fc:a: fe (y)>2 

+ 0--Pk(y))®{(x k (y) + A k (t))log(x k (y) + A k (t))} 

k-x k (y)>2 

A k {t) 1 



J2 Pk(y)^ | (^fe(y) + A k {t) - 1) ( Iog(x fc (y)) + log ( 



1 + 



k-x k (y)>2 



Xk{y) x k (y) 



+ £ (1 -Pfc(l/))E{(x fc (l/) + A fc (*))(log(x fc (y)) + Iog (l + J " (0 



£ (scfc(j/)+E{A fc }-p fc (y))log(x fc (y)) 

fc:z fe (y)>2 



+ £ K fe) E {(M s) + A.( t )-i)io g (i + ^-^)} 



fc:Zfc(y)>2 

+ £ a-Pk(ym\(x k (y) + A k (t))\og[l+ Mt) 

k-x k (y)>2 



Xk{y) 

Now notice that for constants a > 0,6 > 0, c > 



E|(a + i4*(t))log 



= E {(a + A fc (t))(Iog(6c + A fc (*)) - log(c))} 

< E{(a + A k {t))(\og{bc + A k (t))-\og(c))\A k (t) > 1)} 

< E{(a + A k {t))(bc + \og(A k (t))-\og(c))\A k (t) > 1)} 

= a(6c-log(c))+E{6c^(t) + alog(A fc (t))|^ fc (t) > 1)} + E {A k (t) \og{A k (t))\A k (t) > 1} 

< a(bc - log(c)) + i _ a6 p C ^ ( fc t) ( l } 0} + E {i4*(t) Iog(^(t))|A fc (t) > 1} . 
Thus, as E{A fc (t)} and E {A fc (t) log(i4 fc (t))|A fc (i) > 1} are bounded we find 

E|(a + A fc (t))log(6+^)| <oo. (3) 
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Hence AV(y) = G{y) + G 2 (y), with 

G 2 (y) = J2 Pk(y)V{A k (t) > l}E{(l + A k (t)-l)log(l + A k (t)-l)\A k (t) > 1} 

k:x k (y)=l 

+ (l-Pfc(y))E{(l + A fc (t))log(l + A fc (t))} 

k-x k (y)=l 

+ Y, F iMt) > l}E{A k (t)log(A k (t))\A k (t) > 1} 

k-x k (y)=0 

A k (t) 1 



+ E »WE{few + A( t) -i)iog(i + 

+ ^ (l-p fc (y))E{( 3 ; fc (y)+A fc ( i ))log(l + - l ' v(/; 



Xk(y) 



k:x k (y)>2 

which is a bounded function by ^ and as K < oo. 

Lemma 4.3. Assume A/3 < /U. Then, for all states y G -B(e) e = |(1 
G(y)<-en Y v{Pk)\og{x k (y)). 

k-x k (y)>l 

Proof. Since e = |(1 — > we have 

= ^ log(x fc (y)) [(1 - 2e)//i/(i\) - p fc (y) . 

Now note that, for all states y G -B(e), we may write 

Y \og{%k{y))Pk{y) = Y l °s( x k{y)) Y qs ^ y "> 

k-x k (y)>l k:x k (y)>l SeU:S3k 

= Y qs ^ Y l °s(xk(y)) 

S&n k£S:x k (y)>l 

= Y Qs(y)log(ws(y)). 

Sen(y) 



Likewise, we may write 

Y l og(x k (y)) 
k-x k (y)>l 



- Y l °z( x k(y)) Y 1 



k-x k (y)>l 
K 



i:Si(y)3k 



-; X) Y l °s( x k(y)) 

l 1 »_i r..\ 



77 X M^Q/) (?/))• 



□ 



i-Si{y) 
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Substitution of these two equalities in G(y) and remembering that v (P&) = ^ gives 
G(y) = -e/i ^ z.(P fc )log(x fe (y)) 

+ (i-e) X] ^^(^(^(y)) - Y Qs(y)log(ws(y)). 



K 

i-.Si(y) Sen( v ) 



Then, using Lemma 4.1 and recalling the fact that Si(y) G @(y), so that ws i ( y )(y) < w(y) 
for all i = 1, . . . , K, yields 

G(y) < -efi Yl »( p k) l °g( x k(y)) + (l-e)[Y^ lo z( w sdy)(y))- l °g( w (y)) 

k:x k (y)>l i-Si{y) 
k-x k (y)>l 

for all states y G B(e). □ 



By Lemma [42] we know that AV(y) < G(y) + Gf ax , where G™ ax = sup yg ^,G 2 (y) < oo. 
Now consider the set C where the drift of V(y) might be bigger than —1, i.e. consider 

C = {ye*:G(y)>-Gf ax -l}. 



The next lemma shows that this set is small (see Definition A.3). 
Lemma 4.4. Assume A/3 < /i. Then, the set C is small. 
Proof. Consider the sets 

B(e) = (y € * : x k (y) < (S) ' ,Vfc G /C 1 , 



with e = |(1 - ^) > and 

6 = {ye*: logx k (y) < [ 2 + ' ,Vfc G /C 

We see that -B(e) c C 5(e) as subsets of one region are guaranteed and thus w(y) > x k (y), 
for all k = 1, . . . , K . Further we see that C \ B(e) c C C, as follows from the upper bound 
for G(y) found in Lemma 4.3 Hence, 

C C B(e) c U (C\B(e) c ) C B(e) U C = {y G * : x k {y) < M,Wk G /C} , 

where 

M = max e E ^ 



e 

Thus C C Lkm-, where L m is the level set, L m = {y £ * : y(H) < m}. We know by 
Lemma A.3 that Lkm is small and hence, by definition, C is small. □ 
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Using the previous lemmas we can now show stability of our system. 

Theorem 4.5. Assume A/3 < \i. Then, the Markov chain (Y(f)) te fsj is positive Harris 
recurrent with invariant probability measure ir and 



= J n(dx)f(x) < oo, 

where 

fM = ! -G(y)-GT* fory^C, 
J[y> I 1 foryGC, 

Moreover, 

Urn E{g(Y(t))\Y(0) = y} = f n(dx)g(x), Vy G M>, 
/or any function g satisfying \g(x)\ < c(f(x) + 1) for all x and some c < oo. 



Proof. In Lemma A. 3 we have proven that our Markov chain satisfies the irreducibility and 



aperiodicity properties of Theorem A.l Further, we have proven in Lemma |4.4| that the 
set C is small. Thus, as / > 1 by construction and V is nonnegative and finite everywhere, 
we need to show that 

AV(y)<-f(y) + bl c (y), Vy G X I / , (4) 

for some constant b G M, in order to prove our claim. 
For y C we get 



max 



AV(y)<G(y) + G. 

which holds as we have shown in Lemma 14.21 
For y G C we get 

AV{y) <-l + b. 

We therefore take b > 1 + sup y&c AF(y), so that the inequality holds by construction. 
Hence we have shown that Q holds for all y G VP, proving our claim. □ 



Remark. Lemma A. 3 shows our Markov chain to be (/^-irreducible, where ip is the Dirac 



measure on \? assigning unit mass to the empty configuration yo- Hence, by Definition A. 2 



?/>({yo}) > 0. Theorem 3.1 then follows by definition of a Harris recurrent chain, Defini- 
tion \AM 

Remark. The Foster-Lyapunov approach may also be leveraged to derive an upper bound 
for the expected value of functions of the total number of particles in the system. 
Specifically, define J(y) = Efc :£i fe)>i l°g( x k(y)), and denote Gf^ = sup yeB(e)c G(y) + 



J(y) < oo. By virtue of Lemma 4.3 we then have G(y) < — j^J{y). 
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Using Lemma 4.2 and taking expectations yields 

It 
K 



e {v{Y{t + 1))} - e {v(Y(t))} < -i^E {j(y (*))} + G max + cr x , 



for all* = 1,2, .... 

Summing over i = 1, . . . , T, we obtain 



T 

~ X) E { J(F(i))} < E {V(Y(l))} + T(G 1 1 nax + G^ ax ! 
t=i 



and thus 



T->oo J A * — ' 



t=l 



5 Concluding remarks 

We examined the stability of wireless networks whose users are distributed over a compact 
space H. Users arrive at spatially uniform locations with intensity A and each have a 
random number of packets to transmit with mean (5. In each time slot, an admissible 
subset of users is selected uniformly at random to transmit one packet, as governed by the 
prevailing interference constraints. We considered a wide class of interference constraints, 
including the SINR model and the protocol model, and showed that the necessary condition 
A/3 < /i is also sufficient for random admissible-set scheduling to achieve stability, with fi 
denoting the maximum number of users in an admissible subset. 

In the present paper we focused on spaces that are such that any user can always be part 
of a maximum-size admissible set whenever the other users are distributed properly in the 
space. Further we assumed the arrival process to be spatially uniform. Extension to more 
general spaces and non-uniform arrival densities is subject of current research. 
We demonstrated that the necessary condition is not sufficient to achieve stability for seem- 
ingly similar but subtly different scheduling disciplines such as maximal scheduling with 
priorities. An interesting topic for further research is to further demarcate the class of 
scheduling disciplines for which the necessary condition is also sufficient for stability. 
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A Preliminary results 

As mentioned earlier, the stability proof relies on a Foster-Lyapunov approach. In this 
appendix, we first recall various relevant definitions and a result from Meyn and Tweedie [8]. 
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After that, we prove that our Markov chain (Y(t))t£?$ introduced in Section [2] satisfies all 
technical conditions for the Foster-Lyapunov approach to apply. 

Let (Y(i))j e pj be a Markov chain with state space Further, let B(^f) be the cr-field of 
subsets of ^f. This cr-field is assumed to be countably generated, i.e. it is generated by some 
countable class of subsets of iff. 

Definition A.l. (Y(t))i e ^ is said to be (^-irreducible if there exists a measure ip on 
such that, whenever tp(C) > 0, we have 

P |min(t : Y(t) G C) < oo|F(0) = y } > 0, Vy G <3>. 

Let P m {y, C} denote the m-step transition probability to go from state y to the set C G 
B(fy). Further define the transition kernel 

oo 

Ki (y, C) = P m {V, C} 2~( m+1 ), y G 4>, C G B(¥). 

Definition A. 2. (Y(t)) te pj is said to be ^-irreducible if it is y?-irreducible for some <p and 
the measure ip is a maximal irreducibility measure, i.e. it satisfies the following conditions: 

(i) For any other measure cp' the chain is (//-irreducible if and only if ip(C) = implies 

4HC) = o. 

(ii) If ip(C) = 0, then ip{{y : P |min(t : Y(t) G C) < oo|Y(0) = y} > 0}) = 0. 

(iii) The probability measure ip is equivalent to 

ip'(C) = [ <p'(dy)Ki(y,C), 

for any finite measure <p' such that the chain is ^'-irreducible. 

Note. By [8l Thm. 4.0.1] we know that if there exists a measure ip such that the chain is 
99-irreducible, then there exists an (essentially unique) maximal irreducibility measure ip. 

Definition A. 3. A set C G B(^>) is called £ m -small if there exists an m > and a non- 
trivial measure £ m on B(^f), such that 

F m {y,D}>U(D), VyeC,DeB($). 

A set is called small if it is £ m -small for some m > and some non-trivial measure £ m . 

Definition A. 4. Suppose (Y(t))t&i is v>irreducible. The chain is called strongly aperiodic 
when there exists a ^i-small set C with £i(C) > 0. 

Let lc(x) be the indicator function of the set C, i.e. Ic(x) = 1 if x G C and otherwise. 
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Definition A. 5. A ^-irreducible chain (Y(t))t£^ is said to be Harris recurrent if for all 
C G B(4>) such that ip(C) > we have 



Ic(Y(t)) = oo|f(0) =y\=l, VyeC. 



K t=l ) 

If a Harris recurrent chain admits an invariant probability measure it is called positive 
Harris recurrent. 

The following theorem follows from Chapter 14 in [8]. 

Theorem A.l. Suppose that the chain (^(i))teN is ^-irreducible and strongly aperiodic. 
If there exists some small set C , a function f > 1 and some nonnegative function V that is 
finite everywhere such that 

AV(y)<-f(y) + bI d (y), \fy G 4>, (5) 
then (Y(t))t£N is positive Harris recurrent with invariant probability measure ir and 

7r(/) = J ir(dx)f(x) < oo. 
Moreover, 

lim E \g(Y(t))\Y(0) = y\ = f ir(dx)g(x), Vy G 
for any function g satisfying \g{x)\ < c(f(x) + 1) for all x and some c < oo. 



To prove that our Markov chain fulfills the conditions in Theorem A.l, we first show that 
our cr-field, is countably generated. 

Lemma A. 2. is the Borel a-field. Furthermore, B(^>) is countably generated. 

Proof. In this proof we will use some definitions and results in measure theory, see [5] and [6] 
for more details. 

First note that H endowed with the topology generated by the open sets defines a complete 
separable metric space as H is compact. Then it follows by [5| Thm. A2.6.III] that the Borel 
u-field of \I/ is the smallest cr-field with respect to which the map y — > y{B) is measurable 
for any Borel set BOH. Further it follows that ^ endowed with the vague topology is 
a complete separable metric space. The vague topology is the topology generated by the 
mappings (f) —> <pg = J gd(p, with 4> a measure on ^, for all continuous functions g : H — >• M + 
with compact support. 

Since the space is separable, there exists a countable dense set V in this space. Let So be 
the class of all finite intersections of all open sets {x G H : D(x,d) < r}, with d G T> and 
r G Q + . Then, by [5j Lemma A2.1.III], So is countable and generates the Borel a-field, 
B(^). □ 
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We will now prove that our Markov chain satisfies the irreducibility and aperiodicity prop- 
erties of Theorem IA.1I 



Lemma A. 3. Assume A/3 < ji. Then, (Y(t))teN is <f •-irreducible and strongly aperiodic, 
where if is the Dirac measure on ^ assigning unit mass to the empty configuration yo. 
Moreover, the level sets of the form L m = {y G ^ : y(H) < m},m > 0, are small. 

Proof. The proof proceeds along similar lines as in [7]. 

Consider an initial configuration y with y(H) = n < m particles, so y G L m . Remember 
that batches have a size at most one with positive probability. 

First consider the case where batches can be empty. Then P {A = 0} > and the probability 
that the system is empty after m time slots is greater than the probability that no particles 
arrive during the first m time slots times the probability that the n particles are served in 
the first m time slots. Thus, as in a non-empty configuration the probability that at least 
one particle is served in a time slot is at least |, 

F m {y,{y }}>F{A = 0} m Q 

which is greater than zero as ¥{A = 0} > 0. This proves that (Y(t))t^ is (^-irreducible in 
this case, because <f(D) > only when yo S D. 

Now, define the measure £ m = P {A = 0} m (V) m f- For this measure we have P m {y, D} > 
£ m (D) for all D G So we see that L m is £ m -small. Further, {yo} is £i-small and 

£i({yo}) > 0, thus (Y(t))t£® is strongly aperiodic in the case where batches can be empty. 
Now consider the case where batches have a size of at least one and note that the probability 
that a batch has size one is positive. Also note that f3 > 1 and, hence, A < fi. It thus follows 
that P{A = j} > for at least one j G {0, . . . , (J,— 1}. Now denote j* = min{j : P {A = j} > 
0} < [j,— 1, define the initial set of particles by X = {x±, . . . ,x n }, X = IJiLii 2 ' ■ v({ x }) — ^} 
and let Pk, k = 1, . . . ,K be the partition of H given by property [3] of the function F. 
Further let ki be such that Xi G P^ and define hi = [ki/ fi\. Then, define for j = 1, . . . , j*, 
i = 1 n, 



a>j,i-i 



bi/i + j if hoi + j < k 
biH + j + 1 if bill + j > k 



Further, for j = 1, . . . , j* and i > n, denote a^-i = j. 
Now define, 



j 

(t) = J] P {A(t, P a . t ) = l\A(t, P aitt ) = 1, . . . , A(t, P a ._ M ) = 1, A(t, H) = f) 



and note that, for all t, 
K(t) > j^? > 0, 
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as the locations of batches is uniform and independent of the locations of other batches. 
Further note that the probability that a given subset gets removed, given a total of s 
particles in the system is at least ^. Hence, 

P m {y, {yo}} > II IP {A(t -l,H)=j*}F a (t-l) ¥iW >(F{A = f} (2jqj -. 2w j > 0. 

As above, this proves that (Y(t))teN is ^-irreducible, because f{D) > only when yo £ D. 
Now, define the measure MP {A = j*} ^ 2 Ky*2 n ) </3 ' ^ or ^ s measure we nave ^ m {Vi D} ^ 
£m(D) for all D G So we again see that L m is £ m -smail. Further, {yo} is ^i-small 

and £i({yo}) > 0, thus (Y(t)) te ^ is strongly aperiodic. 

□ 
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